Old Dominion University Desearch Foundation 




DEPARTMENT OF MECHANICAL 
COLLEGE OF ENGINEERING 
OLD DOMINION UNIVERSITY 
NORFOLK, VIRGINIA 23508 


ENGINEERING AND MECHANICS 






7S9 



PROJECTION FILTERS FOR MODAL PARAMETER ESTIMATE 
FOR FLEXIBLE STRUCTURES 


By 

Jen-Kuang Huang, Principal Investigator 
and 

Chung-Wen Chen, Research Assistant 


Progress Report 

For the period ended December 31, 1986 


Prepared for the 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 23665 


Under 

Research Grant MAG- 1-655 

Dr. Jer-Nan Juang, Technical Monitor 
SDD-Structural Dynamics Branch 

4BASA-CE-16G3C3) PJBCJ ECTICN EJL1IBS FOB N87-26563 

HCEAL FABAMETEB iSXIMATE FCB E1EX1BLE 

53BUC10BES Progress Eeport, period ending 31 

Eec. 1966 (Old Ecainicn Uri f . ) 30 p Unclas 

Avail: BUS HC AC3/ME A01 CSCL 12A G3/64 0063719 


February 1987 


DEPARTMENT OF MECHANICAL ENGINEERING AND MECHANICS 
COLLEGE OF ENGINEERING 
OLD DOMINION UNIVERSITY 
NORFOLK, VIRGINIA 23508 


PROJECTION FILTERS FOR MODAL PARAMETER ESTIMATE 
FOR FLEXIBLE STRUCTURES 


By 

Jen-Kuang Huang, Principal Investigator 
and 

Chung-Wen Chen, Research Assistant 


Progress Report 

For the period ended December 31, 1986 


Prepared for the 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 23665 


Under 

Research Grant NAG-1-655 

Dr. Jer-Nan Juang, Technical Monitor 
SDD-Structural Dynamics Branch 


Submitted by the 

Old Dominion University Research Foundation 
P. 0. Box 6369 
Norfolk, Virginia 23508 


February 1987 


PROJECTION FILTERS FOR MODAL PARAMETER ESTIMATE 
FOR FLEXIBLE STRUCTURES 

By 

Jen-Kuang Huang 1 and Chung-Wen Chen 2 
ABSTRACT 

Single-mode projection filters are developed for eigensystem parameter 
estimates from both analytical results and test data. Explicit formulations 
of these projection filters are derived using the pseudoinverse matrices of 
the controllability and observability matrices in the general sense. A 
global minimum optimization algorithm is developed to update the filter 
parameters by using interval analysis method. Modal parameters can be 
attracted and updated in the global sense within a specific region by pass- 
ing the experimental data through the projection filters. For illustration 
of this new approach, a numerical example is shown by using a one-dimen- 
sional global optimization algorithm to estimate modal frequencies and 
dampings. 
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Introduction 

The problems of deriving control algorithms and state estimators for 
maneuvering flexible structures have been investigated by many researchers in 
recent years. The control design demands an accurate model of the system 

dynamics which will adequately describe the system's behavior. System 

identification methods use experimental measurements to estimate dynamic 

properties such as natural frequencies, damping factors, mode shapes and modal 
masses which are referred to as modal parameters. Several different time- 
domain and frequency-domain methods are possible for the identification of 
structures. Various techniques may share the same mathematical foundation via 
system realization theory^-. However, most techniques do not account 
explicitly for the factors which will affect the actual performance in 
practice significantly. These factors include nonlinearities, local modes, 
and system and measurement noises. In order to achieve the final purpose of 
identification, i.e., control of flexible structures, an on-line estimation 

technique needs to be developed. This technique may provide updated modal 
parameter estimates only for specific regions needed to be controlled. On the 
other hand, modal parameters can also be identified by using the analytical 

finite element method. The result is usually used only for the comparison 

with the experimental one. However, the analytical result may provide 
valuable initial estimate for modal parameters for an on-line estimator. 

For linear time-invariant systems, optimal model-reduction and state 
estimation has been developed via optimal projection equations based on 

modified Riccati and Lyapunov equations . Other filtering approaches in both 
time and frequency domains are easy to implement and effective in rejecting 
uncorrelated measurement noise from simulated data^. Although filtering 

approaches are not restricted to linear sjystems, time-domain filters usually 
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involve unacceptable computational burden for multi-mode system like large 
flexible structure. 

The objective of this paper is to introduce simple projection filters for 
modal state estimate. These filters are formulated with a single mode only 
and their explicit expressions can be derived using the pseudoinverse of the 
controllability and observability matrices in the general sense. Filter 
parameters are initially implemented from the analytical model and updated by 
real data using a global minimum optimization algorithm. The global minimum 
optimization algorithm is developed by using interval analysis method. Since 
the filters are developed in modal space, system modal parameters within a 
specific region are, as a by-product, identified. 

Finally, a numerical example for a ten-mode structure is given to 
illustrate this new method. A one-dimensional global optimization algoritm is 
also developed and guarantees to find the smallest value of a cost function 
throughout a specific closed region of modal parameters. 

Projection Filters Formulations 

The projection filters are developed based on system realization 
theory. A finite-dimensional, linear, time-invariant dynamic system can be 
represented by the state-variable equations in discrete-time form: 

x(k+l) = A x(k) + B u(k) (1) 

y(k) = C x(k) (2) 

where x is an h-dimensional state vector, u is an m-dimensional control or 
input vector, and y is a p-dimensional measurement or output vector. The 
integer k is the sample • indicator. For flexible structures, the state 
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transition matrix A is a representation of mass, stiffness, and damping 
properties. The control influence matrix B charaterizes the locations and 
type of input control vector u. The measurement influence matrix C describes 
the relationship between the state vector x and the output measurement vector 
y, and characterizes the mode shapes of the system. 

For the state-variable Eqs. (1) and (2) with free pulse response, the 
time domain description is given by the function known as the Markov parameter 

Y(k) = CA k ’ 1 B (3) 

or in the case of initial state response 

Y(k) = CA k x(0) 

where x(0) represents the initial conditions of state vector and k is an 
integer. The functions Y(k) can be obtained from the experimental data and 
used to form the (r+1) by (s+1) block data matrix (generalized Hankel matrix) 


H(k-l) = 


Y(k) 

Y(k+t x ) 

Y(k+t $ ) 


Y( jj+k) 

• 

• 

Y(j^+k+tj) ... 

• 

• 

• Y( jj+k+tg) 

• 

• 

(4) 

Y(J r +k) 

Y(j' r +k+tj) ... 

.Y(j r +k+t s ) 


|(l — 

are arbitrary 

integers. 

For the system 


with initial state response measurements, simply replace H(k-l) by H(k). 
From Eqs. (3) and (4), it can be shown that 


H(k) = V r A* W $ ; V r = 


CA 


CA J r 
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and 


W = [8, A 1 B,..., A S B] 


(5) 


where V p and W $ are generalized observability and controllability matrices. 
The dimensions of V r and W $ are (r+l)p x n and n x m(s+l) respectively. Now 
observe that 

H(0) = V W c (6) 

r s 

we can derive 

V* H ( 0 } W* = I n (7) 

# # 

where Y r and W $ are the pseudoinverse matrices of V r and W $ respectively 

in a general sense. I n is an identity matrix of order n. Now, instead of 

# # 

having the matrix and r o r n-dimensional multi-mode system, we develop 

simpler forms of and which represent the pseudoinverse matrices of 

respective generalized observability and controllability matrices derived from 
a single-mode model only. Note that and W* are rectangular matrices with 
dimensions (r+l)x2 and 2x(s+l), respectively. The general explicit 
expressions of and will be derived later. The matrices and W^, which 
are formulated only for specific modes of interest from the analytical 
results, will be used as the left and right projection filter respectively. 
The Hankel matrix H(0), which is formed by experimental data, will then run 
through the projection filters to attract the system modal parameters. If the 
projection filters have the same modal characteristics as the actual system 
does, then from Eq. (7) we have 

V # H(0) W # = I 2 (8) 

where I 2 is a 2x2 identity matrix. Otherwise, we have 
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V* H ( 0 ) W* = 0 


(9) 


This indicates that the modal parameters of the projection filters are 
different from those of the actual system. The projection filters should be 
tuned in order to match the actual modes. The algorithm for filter update is 
developed in the next section. 


Now, the explicit expressions of the projection filters Y* and.W # can be 
derived as follows. A single mode, continuous-time, linear, time-invariant 
dynamical system has the state-variable equations in modal space 



x = A 

x + 3 

u 

(10) 


<-< 

II 

Ol 

X 


(11) 

wi th 


s 




A = 

a 

-O) 

N 

CO 

a 

> 

(12) 


where u is the modal frequency and a is the damping. The corresponding 
discrete-time system can be represented by Eqs. (1) and (2) with 


and 


A = 



N 


cosujT 

e" aT 

sinwT 

-a T . T 
si nuT 

e -oT 

cosuT 




(13) 


B = 



C = [cj C2 ] 


(14) 


where T is the sampling time and bj , b 2 , c 1# c 2 are scalars. From Eq. (5) 
with j 0 = t 0 = 0 
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V = 


V' s 


CA ' 


[ c 2 + c i ^ (r), - Cj V^(r) + C2 ^2^^] 


t t, t 

W = [A 0 B, A 1 B,..., A S B] = 


> 

-b 2 vj(s) + b 2 vj(s) 

b l V 3 (s) + b 2 V 4 (s) 


where 


T r ~ oT 

V 2 (r) = [..., - e sin(j.uT),...] 


T r ~ j ,- oT 

^2 ( r ) ~ [•••» 6 co s(j\jtoT ),...] 


with i =0, 1 , 2 ,..., r 

T r “ \ <jT 

^3 ( $ ) ~ [•••» “6 s1n(t. wT),... ] 


vies) = [..., e 


-t. aT 


'4' a ' - L • • • » c cos( t^uT) , . . . ] 


with k = 0, 1, 2,...,s 
Assume we choose j 1t t k as follows 

J’r ~ Jr-i = Ji i = 0. 1» 2 ,..., integer [£] 

t s ~ t s-k = ^k k = 0, 1, 2,..., integer [-|] 

Then, the projection filters, V* and W # , may have the following 
expressions (see Appendix for proof): 


(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

expl icit 
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V* = 




r~7 (c 2 v a' (r) * c i < <r>> 
i * c ^ 


C 1 + c 2 


7—7 (_C 1 Y r) + c 2 vj 


(r)) 


C 1 c 2 


(23) 


w# = (-b 2 V*(s) + bj V*(s„ T , -jLj. (bj V*(s) + b 2 V*(s)) T ] (24) 


b l* b 2 


where 


vf(r) 

a 


3^1 sin(j uT) sin ((j -j ) wT ) 

[ -* 6 { ’—T— { T^rJ + ^ ‘ H 


V*(r) 


( ( r ) " IJTrT ) ^ ^ 

j.aT cosCj^T) : , cos((j -j.) u T) 

■[-• * ( 2 ( T^TFT * TJTrT 1 + -T 2 


(25) 


{ X]TrT ~ ^ 2 (r) } } » • • • ] 


(26) 


with i = 0,1,2, . . . ,r 


and 


V C # (S) 


- [. 


t t oT sin(t.uT) . . 

« (- T— (™ * ) - 


TJTsT + TJTsT 


( TJTsT r TJTTF^ ) * * * • J 


sin ((t $ - t k )wT) 


(27) 


#, r t.aT cos( t. ojT) . . 

V <1 (S) '[•••’ e ( 2 ( 50TT + TCTsT 1 * 


cos( (t s -t k )uT) 


' ^3 ( s )" I^TsT 


TTT^ )»•••] 


i 8) 


l 
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with k = 0 , 1 , 2 ,..., s 

X^r) and X 2 (r) are the eigenvalues of V T V, and X-j(s) and X 4 (s) are the 
eigenvalues of ww"*”. X^(r), x 2 (r), ( s ) and X^(s) can be derived as 

follows 


X 1 (r) = 


rl+m+Y ( m) 

*m+Y(m) 


if r is even 
if r is odd 


(29) 


with 


m * { 


r 
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r+1 


X 2 (r) = m-Y(m) 


if r is even 


if r is odd 


m-1 

Y(m) = Z cos((j -2j.)wT) 
i=0 r 1 


X (s) = ( 1+r+z ( n ) 

j s Win) 


X 4 (s) = n-Z(n) 


if s is even 
if s is odd 


(30) 

(31) 

(32) 

(33) 


wi th 


n = { 


s 
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if s is even 


s+1 ... . .. * 

if s is odd 


n-1 

Z(n) = Z cos((t -2t. )ojT) 
k=0 s K 


(34) 


Note that and are rectangular matrices with dimensions 2x(r+l) and 
(s+l)x2 respectively. 

Filter Update 

In order to update the projection filters to attract the actual modes 
from experimental data within a specific range of accuracy, a cost function is 
formed as follows: 


J = U T U 


(35) 
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where 


U - [E u , E 12 , E 21 , E 22 ] 


and 


’ll 

’12 


’12 

'22 


= r h ( o ) w’ - i 2 


From Eq. (8), the cost function J would go to a minimum value (ideally, zero) 
when the projection filters have the same modal characteristics as the actual 
system does. On the other hand, for a specific region of system parameters of 
interest, we may update those system parameters of the projection filters 
within the specific region so that the cost function is globally minimized. 
This global minimum of the cost function in the specific region will provide 
the best estimates for the system parameters of the actual system. Although 
the cost function may be corrupted by the system or the measurement noise, the 
system parameters corresponding to the global minimum are expected to be quite 
insensitive to the noise if the noise is not particularly correlated to those 
parameters we estimate. 


Interval Analysis: Global Minimum Optimization 

The method for computing a global minimum within a specific region of 
system parameters is based on the algorithm developed by Hansen 5-8 . Although 
this algorithm can deal with problems in the multi-variable case with 
inequality constraints 5 , only single variable 5 (either modal frequency or 
modal damping) is considered here. The global optimization algorithm basically 

O 

uses a Newton method 0 in conjunction with the interval analysis to solve a 
system of nonlinear equations. The term "global minimum" used herein refers 
to the smallest valu^ of the cost function J throughout a closed interval of a 
system parameter. Because of the interval analysis, the computatinal 
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procedure of this algorithm requires explicit expressions of the first 
derivative (J') and the second derivative ( 0 ' ' ) of the cost function J shown 
in Eq. (35). This can be easily derived by using the explicit expressions for 
the modal filters shown in Eqs. (23) and (24). The algorithm developed by 
Hansen^ has been slightly modified and summarized as follows: 

Initial Step : The algorithm starts with an initial interval X Q . This 

interval is equally subdivided into subintervals which are stored in a list 
Lq. A list Lj (initially empty) consists of intervals for which the width is 

smaller than a specified value w^ and the corresponding width of J is smaller 

than a specified value Let x denote a feasible approximation to the 

global minimum. If the feasible point is not given, the upper limit of the 
cost function is set to J =® with x indefinite. Let [j^, j R ], 

[j|\ j R ] and [j^ 1 , j R 'J denote the interval resulting from evaluating J, 

J' and J' 1 in interval arithmatic using the argument X, respectively; that is 

J(X) = [j L j R ], J'(X) = [j£, jj], J"(X) - [j'', j''] (36) 

and X = [x|_, x R ] 

Then, use the interval analysis to find the corresponding J, O' and J' 1 for 

all the subintervals in Lq. 

Main Steps : 

1. If the list Lq is empty, go to step 11. Otherwise, find the subinterval 

X in Lq for which the left endpoint of J(X), i.e. j|_, is smallest. 

2. If x e X, set x = x. Otherwise, set x*m(X) = midpoint of X. If 

j|_ > 0, the cost of any point inside the intervals in Lq exceeds the 
upper li it J. Then Lq is set empty and go to step 11. 
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3. Concavity check. 

If ' <0, J is concave in X and cannot have a minimum in the interior 
of X. Then, X is deleted and go to step 1. 

4. Monotonicity Check 

If j R < 0 or Jl > °* the 9 rad ‘ ,en ' t of J i" s strictly positive or 
strictly negative over X. Then, X is deleted and go to step 1. . 

5. Gaussian elimination 

Denote E = 3 - J(x), A = lJ'(x)] 2 + 2Ej^' 

If jV > 0 and A>0, i c implies that J (y ) >3 for any y e X. Then, X 

is deleted and go to step L. Note that this is true only for > 0, 
which is not indicated in ef. 5. 

6. Interval Newton Method^ 

If jV ' > 0, denote S' = x - J* (x)/J‘ ' (X) and S = intersection of S' 

and X. Otherwise, denote S = U S£ 

Here and S2 are defined as follows: 

Denote c * x - J'(x)/j^' for * 0. 

and d = x - J'(x)/j R * for j R ' * 0. 

[x. , d] v/hen jA'>0 and d>x. 

If J'(x)>0 S 1 = { empty when j^., 0 or d<x ^ 

Cc, x R ] when ' <0 and c<x R 

^2 “ f empty when or c>x R 

[x, , c] when j.' '<0 and c>x. 

If J ' ( x ) <0 S x = { empty when Qr c<x ^ 

[d, x R ] when j R ’>0 and d<x R 

^2 ~ ^ empty when ' =0 or d>x R 
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7. If S is empty, go to step 1. If S=X, then split X in half. 

A 

8. For each new generated subinterval X = S 1 or S 2 , repeat steps 3 and 4. 

9. Update 3 

For each new subinterval X, denote 

A A A A A 

w[X] = width of X, x = mid [X] = midpoint of X 

A A A A A A 

X = [x L , x R ] and J ( X) * [j L> j R ] 

If J(x)<3, simply replace 3 by J(x) or conduct a line search to 

reduce J as follows: 

A A A A A 

a. If J ' ( x ) >0 , denote x^ = x^. Otherwise, denote x^ = x R . Set 

A A 

x o = x * 

A A A A A A 

b. Denote x 2 = (Xg + x^)/2. If J ( x^ ) > max [J(Xg), J ( x ^ ) 3 , go to 
step e. 

A A A A A A 

c. If J(Xq) <J(x 1 ), replace x^, by x 2 . Otherwise replace x Q by x 2 . 

d. If I x i” x Q I > IF w ^ x ^» 9° t0 ste P b * 

A A A _ 

e. Set 3 = min[J(x), J(x Q ) , J(x 1 )] and set x to the corresponding 

arugument of 3. 

10. Store new intervals 

A A — A 

For each new interval X, if x^ > J, delete X. 

A A A 

If w[X] < Wj^ and j R - j L < w 2» store X ln L l* Otherwise, store X in 
Lg. Go to step 1. 

11. If the list is empty, go to step 13. Otherwise, delete subintervals X 
for which j L > 3 where J ( X ) = C j R l. 

12. If the list Lj is empty, go to step 13. Otherwise, the midpoints of each 
interval remaining in Lj are used as the global minimums. Note that 
there may exist multiple global minimums. 
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13. Because is empty, the global minimum is located on one of the two 
boundaries of the initial interval X Q , which corresponds to a smaller J. 

Numerical Simulation 

From Eqs. (1), (2), (13) and (14), a linear dynamic system with additive 
measurement noise can be represented by: 


x(k+l) = A x(k) + B u(k) 





y(k) 

>< 

o 

It 

+ v ( k ) 

(37) 

wi th 

✓ 


"N. 





A = 

Al A 
a 2 

• 

• 

• 

0 


B T = [B 

{, B l Bj]. C = [C p C 2 ,.. 

l — l 
& 
O 

• 


0 


A 








m 






l 








-a.T 


-a 

•j 

N. 



1 

i 

e J 

COS U.T 

e 

j 

si n a) .T 



! V 

-a.T 

J 

-a.T 

J 

2 f • • • > rn 

(38) 


-e J 

sin w.T 

e 

j 

cos u.T 





J 



J J 




where w., a. are the modal frequency and damping for mode, m is the 

J J 

number of modes, and v(k) is a white noise. To illustrate applications of the 
projection filters in a single input and single output case, the actuator and 
sensor are chosen and located to give 

B. = [J], C.= [0 1] j = 1,2,. ...m (39) 

From Eqs. (3), (37) and (39) with free pulse response, we have 


Y(k) = CA k " 1 B + v(k) 
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m -a.(k-l)T 

= - 2 e J si n[co .(k-l)T] + v(k) 

j-i J 

This allows one to form a symmetric Hankel matrix H(0) from Eq. (4) by 

r — s , j •j — t'j — i . i - 1,2,. ..,r 
From Eqs. (14) and (39), 

b l = c 2 = 1» t>2 = c 2 = 0 

Then, the projection filters, Eqs. (23)-(28), become 

v # «* - [V.V) T , V>, T ] . 

V* (r) b a 

^ 1 (r) = \ 3 (s), \ 2 (r) * X 4 (s), V*(s) = V*(r) , 
and 

v d (s) = V b (r) 

Because the Hankel matrix is symmetric, the cost function shown in 
can be simplified: 


J 






) 


wi th 


E 11 • V a H(0 > V b T ' 1 “ E 
E 12 * V ! H(0 > V 


(40) 
using: 

(41) 

(42) 

(43) 


Eq. (35) 
(44) 


E 21 ’ K H <°> V b T 
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The first and second derivative of J can be derived as follows: 


j’ = 2 e' u e u + e; 2 e 12 + v 2l 


(45) 


J " * 2 Hi Hi + 2 <Hi> 2 + Hi E i2 + <H 2 > 2 + Hi Hi + <Hi> 2 (46 > 


wi th 

V n * V* H { 0 ) (V^ T ) ' + V* H (0 ) (Vf)' 

E 12 s 2 v a H{0) (V a T) ' 

E 21 s 2 V b H(0) {V b T) ' 

E U a V b H(0) (V a T) " + V a H(0) {V b T) " + 2 (v f > * H (°> (V*V 
E 12 * 2{V a ) ' H{0) {V a T) ' + 2 V a H{0) ^ ' 

Z 21 * 2 (v b } ' H{ °) ( V b T) ' + 2 V b H(0) (V b T) " 

One numerical example is illustrated by using the following conditions 
(See Eqs. (29)-(34)): 

r = s- 49, m * 10, T * 0,04 sec (47) 

w x = 0.001, w 2 = 0.01 

where w^ and w^ are the accuracy of system parameters (either freqency or 
damping) and the accuracy of cost respectively, used in the global minimum 
optimization. Specific modal frequencies with different damping factors and 
noise levels are shown in Table 1. The noise level is the ratio of the noise 
standard deviation with respect to the maximum value of V(k), i.e. the peak of 
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free impulse response. For each case, the simulation starts by forming a 
Hankel matrix for this ten modes structure with a damping factor for all modes 
and a specific noise level. Corresponding to each modal frequency, ten 
frequency intervals are given for the projection filters: 

[.1, 5], [5, 10], [10, 15], [15, 20], [20, 25], [25, 35] 

[35, 40], [40, 50], [50, 60], [60, 70] in rad/sec. 

For each frequency interval, with a fixed zero damping factor, the projection 
filters (Eq. (43)) first update their frequencies by using the interval 
analysis method to find the global minimum of the cost function (shown in Eq. 
(44)) within this frequency interval. The midpoint of the final frequency 
interval (width is smaller than w^) is used for the first estimate of the 
acutal modal frequency. With this estimated frequency, the projection filters 
then update their damping factors with an initial interval from 0 to 12 % by 
using the same interval analysis method to attract the damping factor from the 
Hankel matrix. The midpoint of the final damping factor interval (width is 
also smaller than Wj) is used for the first estimate of the actual damping 
factor. With this new damping factor, the whole procedures are repeated 
again. Since the second estimates of the modal frequency and damping factor 
are quite similar to the first estimates, further estimates are prohibited. 
The percentage errors for the second estimates of the modal frequency and 
damping factor are then calculated for each mode and listed in Table 1 and 2. 
The cost function J is plotted in Fig. 1 as the projection filters update 
their frequencies with a fixed zero damping factor for both cases: no noise 
and 30% noise level with zero damping factor in the data. The result shows 
that the cost function is distorted by the noise, but the minimums are not 

i 

effected too much. From Table 1, as the damping factor varies from 0.3%i to 


17 


10%, the errors of estimated modal frequencies fall within 22 for the noise 
free case. As the noise level increases, the errors increase proportionally 
and stay within 102 for 302 noise. For a fixed noise level, the errors 
increase for most of the modes as the damping factor increases. This may be 
caused by the fact that the signal will damp out faster for higher damping 
factor, especially for high frequency modes. From Table 2, similar results 
are found for the modal damping errors except that the percentage errors are 
higher. For 302 noise, the modal damping errors generally fall within 1002. 
For low frequency and low damping modes, the errors are higher because the 
contribution of the damping is comparably smaller. As a result, this 
numerical simulation shows that the projection filters are a promising way for 
the esimates of the modal frequencies and dampings. 

Concluding Remarks 

Two developments are presented in this paper. First, projection filters 
are developed for modal estimation for dynamical systems. Explicit 
expressions of these single-mode filters are derived using the pseudoinverse 
of the controllability and observability matrices in the general sense. 
Filter parameters are initially implemented from the analytical model and then 
updated to attract the actual system modal parameters within a specific region 
by passing the experimental data through filters. Second, a global minimum 
optimization algorithm is developed by using interval analysis method. Giving 
the first and second derivatives of the cost function, this algorithm 
guarantees to find the smallest value of the cost function throughout a 
specified closed region of system parameters. Numerical simulation with a 
one-dimensional global optimization algorithm shows that the errors of the 
estimates for modal frequencies are less than 22 for low noise case and fall 


18 


within 10% for 30% measurement noise. Several sets of single-mode projection 
filters may be used at the same time for modal parameters estimates in 
different ranges for the control requirements. Multi-variable global 
optimization algorithm needs to be developed to imrpove the estimates for the 
modal parameters. 


Appendix 

For the controllability matrix W and the observability matrix V shown in 

* 

Eqs. (15) and (16), the corresponding pseudoinverse matrices, W* and V*, can 
be derived as follows. First, observe that 


v # V = w W # = I,. 


(Al) 


From Eqs. (17) and (25), it is shown that 

2 

u r sin ( j.uT) , . 

V r> v i (r) * ^ t 2 — (xjmr * ijttt) * 


sin( j i .o)T)sin( ( j r -j i - )«T) 


p 

( ^TrT ' " TX^TrT [ s1 " 2(j l" T) + S1 ' n sfn <U r -j 1 )“ T) ] 


1 T* 2 

+ TZ^rT . Z Q ^ sin (j ’i wT) “ S1 ’ n UjwT) sin ( ( j r ~j i )u>T) ] (A2) 

If r is even, with the aid of Eqs. (21) and (29)-(31), one obtains that 


2 

2 [sin (j.uT) + sin(j.uT) sin( ( j -j . )uT) ] 
i=Q 1 1 r 1 


r/2-1 


2 [sin(j.a)T) + si n ( ( j - jJu>T)]‘ + 2 sin 2 (j coT) 
i =0 1 r i ru 


r/2-1 


= 2 [4 cas*((jj-Np.)oiT) sin 2 -L*-] + 2 sin* (1 uT/2) 

i=0 1 * * r 


j-wT. 


2 
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, r/ 2-1 ? 

2 si n 2 ( j ojT/2 ) [1 + E 2 cos 2 ((j. - J r /2)wT)] 

r i=0 1 r 


9 r/ 2~i 

2 sin 2 (i ujT/ 2) [l + E (1 + cos ((j - 2 j.)cjT))j 

r 1=0 r 1 

2 sin 2 (j toT/2) X^(r) 


E [sin 2 ( j j <*>T ) - s-*'n(j.o)T) sin((j - j,)wT)] 
i=0 1 

r/ 2-1 o 

= E [sin(j.ojT) - si n ( ( j - jJuT)] 
i=0 1 


r/2-1 ? ? 

E [4 cos i j wT/2) sin i ( j . - j /2)coT)] 
i=0 


r/2-1 


= 2 cos 2 ( j uT/2) E [l - co s ( ( j - 2j . )u>T) ] 

I • r\ • * 


= 2 cos^ (j r u)T/2) X 2 (r) 

If r is odd, with the aid of Eqs. (21) and (29)-(31), one arrives that 

E [sin 2 (j^uT) + sin(j.wT) sin((j r - j^JwT)] 
i=0 

( r-l)/2 r l2 

= E [sin(j.u)T) + sin((j - j . ) ojT ) ] 
i=0 1 


(r-l)/2 


= 2 sin 2 (j toT/2) E [l + cos( ( - 2^^. )coT) ] 


= 2 sin 2 ( j p uT/2) X^r) 
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F O - 

Z [sin (j.uT) - sin(j.uT) sin ((j - j.) wT) J 

i=0 1 


(r-l)/2 r . i2 

= Z [sin(j.u)T) - sin ((j - j.)wT)] 

j: _ n * f S 


= 2 cos 2 (j r uT/2) 


(r-l)/2 

2 

i=0 


= 2cos 2 (j r a>T/2) X 2 (r) 


[ 1-cos (( j r -2j 1 )wT) ] 


(A6) 


Substitution of Eqs. (A3)-(A6) into (A2) yields 
V*(r) V 2 (r ) = -jx^TrT ^ 2 sin 2 (j r uT/2) X^r)] + 21 ^ 1 ' rT ^ 2 cos2( J r wT/2)X 2 (r) ] 

= 1 (A7 ) 


Similary 


V*(r) V 2 (r) = 1 

Next, from Eqs. (18) and (25), it is shown that 
u , r 

V*(r) V 2 (r ) = - 2X ~ (r T [ stn ^i uT ^ cosfj^uT) + sin((j r 
, r 

+ TTTrT - 2 sin(j.wT) cos(j.wT) + si n ( ( j p 


(A8) 


- )t»T)cos( j\cjT) ] 


- )coT) cos(j..(dT)] 
(A9) 


If r is even, with the aid of Eqs. (21) and (29)-(31), one obtains 


r 

Z [sin(j.wT) cos(j.uT) + sin((j - j. )wT) cos(j.ojT)] 
1=0 1 
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r/2-1 

E [sin(j.wT) + si n ( ( j - j. )u>T)] [cos(j .u)T) + cos((j - j.)wT)J 
i=0 1 


+ 2 sin (j r wT/2) cos(j r uT/2) 


r/2-1 

Z 4 s i n ( j r ojT / 2 ) cos((j i -j r /2)coT) cos(j r u>T/2) cos( ( j.j-j r /2)u>T) 


r/2-1 2 

sin(j„coT) [l + 1 2 cos ((j. - j /2 ) ojT ) ] 

r 1-0 


r/2-1 

sin(j wT) [l + E (1 + cos((j - 2 j.)wT)] 
r i=0 r 1 


= sin(j r coT) ^(r) 


(A10) 


Z [- sin(jjuT) cos(j-o)T) + sin ((j - j. )ujT)cos( j.wT) ] 

i=0 1 1 

r/2-1 

= Z [sin((j r - j.- )wT) - sin(j.wT)] [cos(j.toT) - cos ( ( j - j. )wT)J 
1=0 1 

r/2-1 

= Z 4 [sin((j r /2 - )(oT)cos( j r uT/2) ] [- si n ( j r ujT/ 2 ) s i n ( ( - J r /2)wT)] 
r/2-1 ? 

= sin(j uT) E 2 sin ((j /2 - j.)wT) 
r 1-0 r 1 


r/2-1 

sin(j wT) E [l - cos((j - 2 j.)uT)] 
r i=0 r 1 
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= sin(j r u>T) \ 2 (r) 

Substitution of Eqs. (A10) and (All) into (A13) yields 
V a (r) V 2 (r) ' [s1n<J r“ T) H (r) ] t s1n( V T) ^ (r) 3 

This is also true if r is odd. 

Similary, it can be proven that 

vJ(r)Y 1 (r) = 0 (A13) 

observation of Eqs. (15), (23), (A7), (A8), (A12) and (A13) leads to 

V* V = I 2 (A14) 

Similar procedures can be used to verify 

WW # = I 2 (A15) 
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Table 1. Percentage Error of Estimated Modal Frequency 
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Table 2. Percentage Error of Estimated Modal Damping 
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Figure Legend 


1 -Numerical example for the cost function J as a function of modal 
frequency with a zero damping and two different measurement noises 
(0% and 201 ) for a ten-mode structure. 
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